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Abstract 

This paper presents an electromagnetic velocity (EMV) 
gauge technique which allowsuse of conductive samples in 
the measurement of material velocity subject to shockwave 
compression. A dynamic model of the induced current in the 
EMV gauge caused by the conductive samples moving in a 
uniform magnetic field is developed and the motional emf 
generated in the EMV gaugeis extracted from the measured 
signal. Experimental results are presented and discussed to 
demonstrate the effectiveness of the technique. 
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Introduction 

Piezo-resistive film gauges are used regularly to 
determine normal stress in matter under shock 
compression. The functionality of these gauges relies 
on a pre-calibrated relationship between the change in 
relative resistance and the change in normal stress. 
Under extreme conditions, both substrate and 
adhesive constituents of these gauges can undergo 
phase and/or chemical changes whose effects on the 
overall performance of the gauge are not clear 
(Rosenberg, 2009). Electromagnetic velocity (EMV) 
gauges are another class of film gauges which permit 
the direct in-situ measurement of shocked material 
flow velocity. The active sensing element does not 
need pre-calibration; rather, this method requires 
exposure of the moving element to a well quantified 
external static magnetic field in order to produce 
motional electromotive force (emf). While applications 
of this technique for measuring velocities of shocked 
materials have been widely reported in the literature 
(Sheffield, 1999; Fritz, 1973),it is limited to non- 
metallic materials only. This is due to signal distortion 
caused by the induced currents within electrically 
conductive materials moving in an externally applied 
magnetic field at high velocities. To use an EMV gauge 


with conductive samples, the motional emf generated 
in the pickup foil has to be extracted from the 
measured signal which results from the superposition 
of both motional emf and voltage from induced 
currents. 

In this paper an electromagnetic technique is 
developed which uses an analytical model to estimate 
the dynamics of induced current between the copper 
substrate and the EMV gauge. The former is modelled 
as a non-magnetic conductor moving in a uniform 
magnetic field with a uni-directional translational 
velocity. The Foucault current in the moving 
conductor is discretized as a piece-wise constant signal. 
For a mm thick copper sample disk moving in a static 
magnetic field with a uniform translational velocity, 
the volume conductor is modelled as a magnetic 
dipole loop carrying a current and the EMV gauge is 
modelled as a closed circuit loop moving in the field of 
the magnetic dipole. The field produced by the 
magnetic dipole is calculated throughout the space in 
which the pickup foil gauge moves. Equations of 
mutual induction are derived and the induced current 
in the EMV gauge loop is solved numerically, which 
allows separation of the emf signal for measurement 
of the particle velocities. Numerical analysis is 
provided for the induced current in the EMV gauge 
with respect to the Foucault current in a copper 
sample undergoing a step wise change in velocity. A 
few experiments were performed in an attempt to 
demonstrate the feasibility of the proposed technique 
to compensate for signal distortion in measuring 
particle velocities of shock compressed copper. The 
ability to use EMV gauges within samples of high 
conductivity complements our use of manganin film 
stress gauges to quantify shocked states in matter. 

EMV Gauge Modelling 

Fig. 1 shows the setup for the EMV gauging technique 
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described in this paper. A shock wave is transmitted 
into a sample, comprising a copper disk above and 
below an insulated metallic foil element, where both 
copper and foil element are set in motion and, under 
the influence of an external magnetic field, generate a 
signal that is proportional to their velocity. According 
to Faraday's law, the particle velocity of the sample is 
given byi; = Bl/V where V 1 is the voltage across the 
active foil element, B is the average magnetic flux 
density as measured with a Gaussmeter at the location 
of the gauge, l is the length of the active element and v 
is the velocity of the gauge. 


r = (a 2 + y 2 + z 2 — 2ay sin <p) 1/2 (3) 

The differential current loop element satisfies 


d\ = —a sin cp dcpi + a cos cpdcp] (4) 

The cross product d\ x rbecomes 


d\ x r 


J 


a sin cpdcp a cos cpdcp 0 
—acoscp y — asincp z 
adcp [z cos cp i + z sin cp j + (a 


(5) 


■ y sin cp) k] 


From eq. (1) the magnetic field at a point P on the 
EMV gauge loop surface is 


Shock generated via explosives or flyer impact 



FIG. 1 EMBEDDED ELECTROMAGNETIC VELOCITY GAUGES 
WITH EXTERNALLY APPLIED MAGNETIC FIELD 


A Model of the Induced Magnetic Field 

In this paper, only the copper disk above the gauge is 
modeled as a loop of radius a carrying a current Ii and 
moving along the z-axis with a constant translational 
velocity v toward the EMV gauge, as shown in Fig 2. 
The effect of induced current from the copper disk 
downstream from the gauge is for now neglected for 
simplicity. 



By the Biot-Savart law, the magnetic field at a point P 
on the surface of the EMV gauge loop is 


where 



d\ x r 

y 3 


(1) 


r = — acoscpi + (y — asincp )j + zk (2) 

with a magnitude 


B x = 0 (6) 

Mo /i az f 2n sin cpdcp 

An J 0 ( a 2 + y 2 + z 2 - 2ya sin <p) 3 / 2 ^ 

= ha f 2n (a -y sin cp)dcp 

z An J 0 (a 2 + y 2 + z 2 - 2yasin<p) 3 / 2 

B y and B z can be expressed as elliptic integrals and 
computed numerically. On the z-axis where y=0, B y 
and B z are 


By = 0, B z 


iM lQ 2 
2(a 2 + z 2 )2 


(9) 


Models of the Copper Disk and EMV Gauge Loops 

Let the mutual inductance of the copper disk and the 
EMV gauge loopsbeM 12 . The circuit equations for the 
two loops can be expressed as 


dB dM 12 dl 2 

RlIl + Ll lF + ^ir l2 + Ml2 ^F = Vl 


dl 2 dMi? dl i 

^2 -rf + R2I2 + h + M 12 -4 = 0 


dt 


dt 


dt 


(10) 

( 11 ) 


In a configuration where the distance between the 
copper disk and the EMV gauge is much smaller than 
the radius of the current loops (i.e., z « a ), the 
magnetic field generated by the source loop can be 
approximated by the magnetic density on the z-axis, 
given by eq. (9), and the magnetic flux going through 
the EMV gauge loop can be simplified as 


. 9. Mo nfaf'b 2 

M n h = B z (0,0,z)(nb 2 ) = 3 

2 (a 2 + z 2 )2 

And the mutual inductance can be expressed as 
Ho na 2 b 2 

~ 3 

2(a 2 + z 2 )2 


(12) 


(13) 


Thus the total time derivative of M 12 is 
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dM 12 dM 12 3p 0 na 2 b 2 v z 

~dT = ~d^~' V= 2 (z 2 + a 2 ) 2 - 5 (14) 

wherep is the moving velocity of the copper disk 
toward the EMV gauge. To simplify the problem, we 
assume that the current running around the copper 
disk loop b holds a constant value in a short time 
interval so that eq. (11) becomes 


dl 2 dMi 2 

L 2 -r + R 2 h = — -r^/i 

z dt z 2 dt 1 

Applying eq. (14), we get 

dl 7 3p 0 na 2 b 2 v 


L 2 + Rih — h 


(15) 


(16) 


2 (z 2 + a 2 ) 2 - 5 

To perform a modal analysis for the induced current I 2 , 
we define the following dimensionless ratios 


v L 2 v I 2 2 aL 2 

t = -t f p = — , / = A-,A = (17) 

a aR 2 b 3p 0 nb z 

Substituting (17) into (16), we derive the modal 
equation for the induced current in terms of the ratio / 


dl 1 r 

dr p (t 2 + l) 2 - 5 v 

which can be solved numerically as shown in Fig. 
3. From the numerical solutions of the induced current, 
we have observed that 


a. The Foucault currents in the copper disk moving 
with a uniform velocity induce currents on the 
EMV gauge before the gauge starts to move. 

b. The induced current causes a reversal of the emf 
signal with a positive overshoot. 

c. The positive overshoot extends for a longer 
timewith higher copper disk velocity. 



Time Ratio (t) 

FIG.3 MODES OF INDUCED CURRENTS IN THE EMV GAUGE 


Modal Analysis 

The procedures for applying modal analysis to the 


experimental datafor signal extraction are as follows: 

• Stepl: Use a Fast Fourier Transform to low-pass 
filter the raw experimental data in order toget a 
smoothed curve so that the mode shapes are 
evident. 

• Step2: Determine the boundary values of the emf 
that corresponds to the modes of the induced 
current demonstrated in Fig. 3. 

• Step 3: Solve eq. (18) as a boundary value problem 
and obtain a modal solution. 

• Step 4: Apply the modal to the smoothed 
experimental data by means of superposition. 

• Step 5: Repeat Steps 2 to 4 until the mode shape is 
not evident on the curve of the experimental data. 

• Step 6: Combines all the modals solved from Steps 
2 to 5 to construct a complete induced emf profile 
and extract the motional emf from the raw signal. 

Experimental Tests, Results and Discussion 

In order to test the performance of the models of 
signal extraction for EMV gauging in shock wave 
environments, a number of trials were performed 
whereby composition C-4 military explosive, hand 
packed in a PVC tube (52mm i.d. x 4.2mm wall x 
203mm high for test 275A and 20.2mm i.d. x 3.3mm 
wall x 101mm high for test 295 A) was used to generate 
an incident shock wave in copper. 

The EMV gauges used for the tests were manufactured 
in-house (Fig. 4 (a)). A 25pm thick brass foil was cut 
with an active length of 10 mm, which produced 
motional emf when moving through a magnetic field. 
The foils were encapsulated with 250 pm (test 295 A) or 
500 pm (test 275 A) thick Teflon adhesive sheets on 
both sides. The completed gauges were then 
embedded between two copper disks (31.75 mm dia. x 
3.09 mm thick above and 22 mm thick below for test 
295A, and 57.05 mm dia. x 3.09 mm thick above and 10 
mm thick below for test 275A) using liquid epoxy 
adhesive, as shown in Figs. 4(b) and (c). 


FIG. 4 (A) IN-HOUSE MADE EMV GAUGE AND (B) EMV GAUGE 
IN THE TARGET SAMPLE (C) DONOR/TARGET ASSEMBLY 
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For trial 275A, a 75 Q resistor was connected in series 
with the EMV gauge and terminated with 75 Q on a 
single scope channel for impedance matching of the 
whole gauge circuit loop. For trial 295A, each lead of 
the EMV gauge was connected to a different 
oscilloscope channel via the center conductor of a 
1.5 m long RG-179 coaxial cable, resulting in a 
differential input to reject common mode noise. Each 
channel was terminated with 13 and the voltage 
difference between these channels provided us with a 
measure of voltage across the foil, accounting for 
resistive losses across cables. The top and bottom 
copper disks were grounded to a central grounding 
point at the instrumentation bunker for all tests. 

An electromagnet built in-house was used to provide 
an external static magnetic field. The magnet 
essentially consisted of 50.8 mm square A36 mild steel 
bars used as pole pieces, two of which were wrapped 
with #8 AWG cable. A switching DC power supply 
was used to provide the coils with a constant current. 

Calibration of the electromagnet was performed before 
every trial. The variation of magnetic flux density (B) 
with vertical position along the center of the pole piece 
gap where EMV gauges werelocated was measured, 
and an average flux density was used for simplicity. 
The variation in B along the sensing width of the EMV 
foil was neglected, since lateral fluctuations in the 
magnetic field were within better than 1% over the 
gauge width. 

The EMV gauge model developed in this paper is 
applied to the velocity data of the copper samples 
acquired in two of the EMV gauging trials. A separate 
experiment was performed to determine the normal 
stress at the location of the EMV gauge within the 
target sample for trial 295A using a Dynasen Inc. 
MN10-0.050-EPTFE manganin gauge. The copper 
material velocity was then calculated via the known 
Hugoniot data (Marsh, 1980)given the normal stress at 
the EMV gauge location. The measured peak normal 
stress was 20.6GPa which corresponds to a particle 
velocity of 0.50mm/psec. 

From the raw experimental data (Figs. 5 and 8), we can 
observe that in both trials the waveforms first show a 
negative peak with a magnitude ranging between 
0.34-0.41mm/psec. The EMV gauge model developed 
in this paper suggests that the moving top copper disk 
perturbed the magnetic field and induced a current in 
the EMV gauge circuit before the gauge produced 
motional emf. This is followed by a positive peak. 


wherethe motional emffrom the EMV gauge is 
embedded. 

The modal analysis is performed on both waveforms 
using the procedures described in last section. 
Twenty-six modals (Fig. 6) for the data of trial 295A, 
and thirty-five modals (Fig. 9) for trial 275A, were 
solved at various time instants with suitable strength 
(i.e., satisfying the boundary values), and then were 
applied to the waveforms to compensate for the 
negative signal distortion and positive signal 
perturbation on the EMV gauge. It should be noted 
that all the EMV gauge model outputs are shown in 
terms of velocities that are converted from the 
corresponding induced current in the EMV gauge. The 
results are shown in Figs. 7 and 10. The negative 
reversal signals were nearly eliminated and the 
positive signals were also compensated for the 
extended effect of the induced currents. In Fig. 7, the 
extracted EMV gauge signal shows an average peak 
value of 0.53mm/psec for trial 295 A, which matches 
the peak value of 0.50 mm/psec inferred by the 
manganin gauge for trial 295 A. The charge diameter 
for trial 275A was about 2.6 times larger than and 
twice as long as the donor charges in all other tests. A 
higher average peak value of 0.72mm/psec particle 
velocity is observed from the output of the EMV gauge 
model for trial 275 A, as shown in Fig. 10, which 
appears to confirm the expected greater delay in both 
lateral and rear expansion waves. These values 
suggest that the EMV gauge model provides a correct 
measure of the incident particle velocity. 

It can be also observed in Figs. 7 and 10 that the EMV 
gauge model increases the period over which an 
effective measurement of particle velocity can be made. 
In trial 295A the effective measurement time was 
extended from about 0.15 to 0.26 psec, whereas in trial 
275A it was extended from about 0.2 to 1 psec. 


0.6 
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FIG. 5 PARTICLE VELOCITYRAW DATA CONVERTED 
FROMEMFMEASURED BY THE EMV GAUGE IN TRIAL 295A 
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FIG. 6 CONVERTED MODES OF INDUCED EMF SOLVED FROM 
THE INDUCED CURRENT MODEL FOR TRIAL 295A 
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FIG.7 COMPARISON OF THE INPUT, INDUCED AND 
EXTRACTED EMV GAUGE SIGNALS FOR TRIAL 295A 
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FIG.8 RAW DATA OF THE PARTICLE VELOCITY CONVERTED 
FROMEMF MEASURED BY THE EMV GAUGE IN TRIAL 275A 

0.12 r 
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FIG.9 CONVERTED MODES OF INDUCED EMF SOLVED FROM 
THE INDUCED CURRENT MODEL FOR TRIAL 275A 



FIG. 10 COMPARISON OF THE INPUT, INDUCED AND 
EXTRACTED EMV GAUGE SIGNALS FOR TRIAL 275A 


Conclusions 

A model was developed in attempt to discriminate 
motional emf from the resultant signal produced by an 
electromagnetic velocity (EMV) gauge embedded 
within shock compressed copper. While the current 
model only accounts for the effects of bulk motion of a 
copper sample on the upstream(relative to the 
direction of shock propagation) side of the EMV gauge, 
the results suggest that it is possible to eliminate the 
negative signal distortion induced by a moving 
conductor in an external magnetic field, thereby 
extending the measurement of particle velocity history 
(up to 1 psec for the cases examined in this 
paper). Although the EMV gauge technique developed 
in this paper requires further development to account 
for the fact that a moving conductor was also located 
dowstream of the gauge, we expect that a similar 
model designed to capture a bulk copper volume 
moving away from the foil gauge could further extend 
the useful part of the velocity signal. Further model 
refinement will be needed to confirm this. 
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